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._!. ■ Abstract. We present preliminary results on the resummation of leading and next- 

T-H , to-leading logarithms for the thrust distribution in deep inelastic scattering. Our 

predictions, expanded to ©(ag), are compared to corresponding results from the 
^ '. Monte Carlo programs DISASTER++ and DISENT. 
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2 ■ PACS numbers: 12.38.Bx, 12.38.Cy 
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0\ ■ 1. Introduction 

C^! For some time now it has been standard practice in e~^e reactions to compare event- 

Qh! shape distributions with resummed perturbative predictions (see for instance [|^, 0). 

r-| ! The resummation is necessary because in the two-jet hmit (small values of the shape 

>■ ! variable) the presence of large soft and coUinear logarithms spoils the convergence of 

k> , the fixed-order calculations. Such resummed analyses have led to valuable information 

about the strong coupling constant and also about non-perturbative effects. 

At HERA, studies of event-shape distributions are being carried out by both 
collaborations but as yet no perturbative resummed calculations exist for comparison. 
Here we present preliminary results on such a calculation Q]. 

In DIS, event shapes are defined in the current hemisphere of the Breit frame to 
reduce contamination from remnant fragmentation which is beyond perturbative control. 
The distribution of emissions in the current hemisphere (Tic) is analogous to that in a 
single hemisphere of e~^e~ . Differences arise from e'^e' however because the momentum 
of the current quark (the quark struck by the photon) depends on remnant-hemisphere 
{'H-r) emissions through recoil effects. This necessitates a resummed treatment of 
the space-like branching of the incoming parton. Currently such a simultaneous 
resummation of the space-like and time-like (double logarithmic) contributions exists 
only for jet multiplicities |0, or for cross sections at large x [0. 

t Talk presented by V. Antonelli at the UK Phenomenology Workshop on Collider Physics, September 
1999, St. John's College, Durham. 
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For the thrust T (as defined in the following section) the recoil from remnant 
hemisphere emissions can be divided into two parts: a piece from soft and collinear 
emissions, which gives double logarithms a" In "^(1 — T), identical to those from half 
a hemisphere of e~^e~, and a piece from purely collinear emissions which gives single 
logarithms a"ln"'(l — T) that can be identified with a change of scale in the parton 
distribution. This is outlined in the next section. 

A valuable cross-check of any resummation is to expand it to NLO and compare 
the result to that from fixed order Monte Carlo programs such as DISASTER++ [§] 
and DISENT |Q. This is illustrated in section |. 

Finally in section §, we comment on possible future developments. 

2. The Thrust in DIS 

There are several possible definitions of the thrust in DIS, which differ according to a 
choice of axis and the normalisation. We here consider the thrust with respect to the 
photon axis, fi^, and normalised to Q/2, 

T=^5^fc:-n,, (1) 

where the sum extends over all particles in the current hemisphere. He- At lowest order 
T = 1, so the region requiring resummation will be that of 1 — T close to zero. 
Expressing all momenta in terms of their Sudakov components 

ki = aiP + /3iP' + kt,i , ai/3i = kyQ^ , (2) 

where P = |Q(1, 0, 0, — 1) and P' = |Q(1, 0,0,1) are along the remnant and current 
directions respectively, we can write the thrust as |]TD[ 

1 - T = ^ A + 5^ a, + a, , (3) 

where the sums now run over emitted particles only and a^ is the a component of the 
current quark^. Dependence on the remnant hemisphere emissions arises both through 
the ^jgT^ Pi term and through 






J2kt, 



(4) 



(we have made use of the fact that f3q ~ 1). 

We now consider the cross section cr(r) for 1 — T to be smaller than r. Roughly it 
will contain virtual corrections from the exclusion of all emissions whose contribution 
to (^ is larger than r. 

§ The expression must be modified when the current quark goes into the remnant hemisphere — but 
such a situation is not relevant for smaU 1 — T. 
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In the remnant hemisphere we exclude soft and colhnear emissions, kt > (3 > t, 
where the kt > (3 condition ensures that particle be in Ti^. The integration over kt and 
j3 gives a double logarithm. In He we exclude the soft and collinear region kt > a > t. 
Again this gives us a double logarithm. The exclusion a„ ~ A;^^ > r is just a limit on 



the maximum emitted transverse momentum: this implies 'stopping' the DGLAP [|Tl 
parton evolution at a scale tQ^, which gives single collinear logarithms. 
These aspects are illustrated in the first order result for a{r): 



a ~ 2tt 



2 In^ - - 3 In - 

r r 



z) 
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(f 



<l(-]Pmi^)+9(-]P,,{z) . (5) 



Os Inl/r f^ dz 
2vr q{x) J^ z 

The first line contains the soft and collinear double logarithms (which turn out to 
be identical to the e~^e~ result) from the conditions on the a^ and /?j. The second 
line contains the collinear logarithm associated with the restriction on the DGLAP 
evolution. Photon-gluon fusion plays a role only through this collinear logarithm, in the 
convolution with the gluon distribution. 

The actual details of the resummation will be presented in Q] . Schematically, the 
result is 



^ = [l + «sC(x)]^^^2^S(r), S(r)=e-^[^''^^^-^''^^M"^^^^^), (6) 



where S(r) is just the corresponding e~^e~ quantity of |I| and C{x) is a r-independent 
coefficient function which is given in [Q. The ratio (^q2\ comes from the suppression 
of radiation with k^ > tQ"^, as mentioned above. 



3. Comparison with fixed order programs. 

Expanding (|]) to NLO it is possible to perform a comparison to the fixed order Monte 
Carlo programs DISASTER++ g] and DISENT 0. We actually look at the coefficient 
of (as/2vr)2 in 



r da{T) 
ao dr 



(7) 



where we have normalised to the the Born cross section ctq for simplicity. At second 
order eq. (|^) contains terms In" r, n < 3 and for the resummation to be correct to 
next-to- leading order we should correctly control terms 1 < n < 3. So the difference 
between our expanded result and the exact result should at most be a constant for 
small r. The upper two plots of figure |I| show our results for the quark and gluon- 
initiated components of the answer compared to the predictions from DISASTER++. 
The shape of the distributions is well reproduced for small r. The lower two plots show 
the difference between the DISASTER++ results and ours: one sees that for small r 
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this difference is indeed compatible with a constant, as required. Comparisons have also 
been made with DISENT, which seems to disagree with our result in the gluon sector 
at the level of a term proportional to In r, and perhaps also in the quark sector. 
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Figure 1. Coefficient of (#^) in -^^. See text for details 



4. Outlook 



Although we have results for the resummation of the thrust distribution in DIS, some 
work remains to be done for the practical implementation of our results for comparison 
with experimental data. In particular, prescriptions need to be defined for the matching 
of our resummed result with the fixed order results, in order extend the range of 
applicability of the calculation (without the matching the results are applicable only 
to very small r). 
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Work is also in progress on the resummation of other DIS event-shape variables. 
Among these, the thrust normalised to the energy in the current hemisphere is close to 
completion. Other variables to be studied include the jet mass, the C-parameter and 
the jet-broadening (the resummation of the latter is relevant also for predicting the form 
of the power correction [Q). 

Once this programme is complete, we hope that it will lead to phenomenological 
studies analogous to those already being carried out in e~^e~ . 
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